Fabrication and Characterization of Topological Insulator Bi 2 Se 3 Nanocrystals 

S.Y.F. Zhao 1 , C. Beekman 1 , L.J. Sandilands 1 , J.E.J. Bashucky 1 , D. Kwok 2 , N. Lee 2 , A.D. LaForge 3 ,S.W. 
Cheong 2 and K.S. Burch 1 B 

1 Department of Physics & Institute of Optical Sciences, University of Toronto, 60 St. George Street, Toronto, 
ONM5S1A7 

2 Rutgers Center for Emergent Materials and Department of Physics and Astronomy, Rutgers University, 
136 Frelinghuysen Road, Piscataway NJ 08854, USA. 

3 Department of Physics, University of California, Santa Cruz, 1156 High Street, Santa Cruz, California 95064, 
USA. 

(Dated: 11 April 2011) 

In the recently discovered class of materials known as topological insulators, the presence of strong spin-orbit coupling 
causes certain topological invariants in the bulk to differ from their values in vacuum. The sudden change of invariants 
at the interface results in metallic, time reversal invariant surface states whose properties are useful for applications 
in spintronics and quantum computation. However, a key challenge is to fabricate these materials on the nanoscale 
appropriate for devices and probing the surface. To this end we have produced 2 nm thick nanocrystals of the 
topological insulator Bi2Se3 via mechanical exfoliation. For crystals thinner than 10 nm we observe the emergence 
of an additional mode in the Raman spectrum. The emergent mode intensity together with the other results presented 
here provide a recipe for production and thickness characterization of Bi2Se3 nanocrystals. 



Topological metallic surface states are predicted to have nu- 
merous properties that are useful for spintronics and quantum 
computation. 1*2 A challenging aspect of this research has been 
to isolate surface state contributions to the measured proper- 
ties of topological insulators 3-7 . Studying nanometer thick 
crystals allows one to tune the chemical potential with elec- 
tric field 8 -^ 2 -, as well as observe modifications of the exci- 
tation spectrum produced by interactions of top and bottom 
surface o 11 ' 13 ^ . One option for producing thin crystals is me- 
chanical exfoliation. This method, together with Raman spec- 
troscopy, has proven to be extremely fruitful in the study of 
graphen o 17 ' 18 and other nanocrystals 19-20 . Indeed, Raman pro- 
vides direct access to the phonon spectrum and can be used to 
map the thicknes s 17,21 or doping level 22 over a large area. 

Similar studies of the topological insulators Bi2Se3 and 
Bi2Te3 have been attempted. To date, these experiments have 
been limited to exfoliated Bi2Se3 crystals > 10 nm thic k 12 ' 23 ' 24 
or Bi2Te3 where the bulk gap is small 8 ' 9 . A limiting factor in 
these experiments was the strong optical absorption of these 
compounds, making the identification of thin crystals difficult 
on Si/SiC>2 substrates. With this in mind, we have mechani- 
cally exfoliated Bi2Se3 crystals on Mica, enabling us to op- 
tically identify crystals only 2 nm thick. By systematically 
studying these crystals with Raman spectroscopy and opti- 
cal and atomic force microscopies (AFM) we have devised 
a method for characterizing the thickness of Bi2Se3 nanocrys- 
tals using non-invasive, all-optical methods. Specifically, an 
additional mode appears in the Raman spectra for ultrathin (< 
10 nm) crystals. The observed thickness dependence of the 
emergent mode intensity can be used for thickness verifica- 
tion of nanocrystals via Raman measurements. 

Bi2Se3 forms a rhombohedral lattice in which the unit 
cell is composed of three five-layer stacks known as quintu- 
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pie layers (QL). A unit cell measures 2.87 nm along the c 
axis, and 1 QL is w 0.96 nm thick. 25 Atoms are arranged 
into planar hexagonal sheets with the sequence -[Se( 2 )-Bi- 
SeW-Bi-Se^ 2 '']— . The superscripts indicate the structural in- 
equivalence of the Se ions, with the SeW atom at a center of 
inversion within a QL. Therefore one expects phonon modes 
to be exclusively infrared (IR) or Raman active. The weak 
van der Waals bonds between neighboring Se^ 2 * 1 planes en- 
ables mechanical exfoliation. Ultrathin Bi2Se3 nanocrystals 
could be identified after exfoliation by optical microscopy in 
transmission mode. The crystal thickness was subsequently 
determined using AFM 26 . 

Raman spectroscopy (A = 532 nm, spot size ~ 1 /im) con- 
firmed that these nanocrystals are indeed Bi2Se3 and was used 
to study the evolution of the crystal lattice structure with thick- 
ness. A raman spectrum (I(uj)) for a bulk crystal is shown in 
Fig- Qji, where we observe a strong phonon mode at 175 cm -1 
and the onset of a mode below 150 cm -1 . Previous results 25 
suggest the higher energy mode corresponds to a Raman ac- 
tive Ai g mode and the latter to an E g mode, in accord with 
group theory predictions for phonons at the Brillouin zone 
center (q = 0) probed by optical experiments. Indeed, the 5 
atoms in the unit cell should lead to 12 optical modes, each 
of which is either exclusively Raman or infrared active^ 5 .. We 
also observe a broad shoulder between 200-350 cm -1 , simi- 
lar to a feature observed in a recent IR study of bulk Bi2Se3 
(see FigEh).^ The presence of this feature in both Raman as 
well as IR and its broad lineshape suggest that it is due to two- 
phonon excitations. 

To quantitatively analyze the evolution of the spectra, we fit 
the measured Raman data with multiple Lorentzian oscillators 
in the form: J(w) = 7 + E i ( (i^B^+f^y ) where i ranges 
from 1 to 3 or 4, depending on thickness, Iq accounts for the 
background, Ei is the center, Ti is the width, and Ai is the area 
of peak i. The resulting fit for the bulk spectrum with three os- 
cillators is shown in Fig. Qji, where the center frequency of the 
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FIG. 1. (Color online) Raman spectra for Bi2Se3 (a) Spectrum for the bulk 
crystal (line + symbols) with the corresponding fit (solid line) consisting of 
three Lorentzian oscillators (dotted lines), (b) Spectra for crystals of varying 
thicknesses, 3 — ► 16 QL and bulk, (c) Spectrum for 3QL nanocrystal. The 
data (line + symbols) is fit with four Lorentzian oscillators (solid line). The 
fit resulting from three oscillators is shown for comparison (dashed dotted 
line). 



E g mode (Peak 1) was fixed at 131.5 cm -1 based on previous 
studies 25 . Turning to Fig. QJ> we examine the evolution of 
the Raman spectra with varying crystal thicknesses, where an 
enhancement of the overall signal is observed with decreasing 
number of QL. As is discussed below this enhancement results 
from multiple reflections in the Bi2Se3 crystal. Perhaps more 
surprisingly, an additional mode emerges at 158 cm -1 as the 
crystal thickness is reduced (see Fig. [TJ:). Indeed, the use of 
only three oscillators, which worked well for crystals thicker 
than 10 QL (Fig. [T^), results in a large difference around 158 
cm -1 between the fit (dashed dotted line) and the spectrum. 
By simply adding another mode (i = 4), the fit (solid line) 
agrees very well with the data. The emergence of the mode 
at 158 cm -1 for ultrathin crystals suggests it can be used to 
verify sample thickness. 

We explore this possibility in Fig. [2] where we plot the 
ratio of the emergent mode intensity with the main peak in- 
tensity (4^). For thicknesses below 10 QL this ratio increases 
and agrees well with ^q^_ 2 ) behavior (dashed line), which 
corresponds to the relative weight of surface and bulk modes. 
Indeed, the intensity of a particular mode is proportional to 
the volume over which the light can emanate. For a surface 
mode, this volume is independent of thickness. In contrast, 
for a bulk mode the volume probed is proportional to the crys- 
tal thickness (for thicknesses less than the penetration depth). 
Unfortunately, due to the detection limit of our instrument, 
we are not able to convincingly detect the mode for A4/A2 < 
0.02. 

While the mode appears to have its origin in the surface, 
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FIG. 2. (Color online) (a) Ratio of emergent peak intensity to the main 
peak intensity (determined from fits) as function of crystal thickness. The 
dashed line is a guide to the eye and plotted to show the ^q^_ 2 ) tren d °f 
the ratio. The gray area for A4/A2 <0.02 indicates the detection limit of our 
system. The ratio A4/A2 as function of crystal thickness due to FP effects 
(solid line, right axis), (b) Left axis: Raman spectrum for a 3QL thick crystal 
(line + symbols) and the lorentzian corresponding to the emergent mode (dot- 
ted line). Right axis: optical conductivity data (ssolid line) of a bulk Bi2Se3 
single crystal taken at T = 5 K. 



it could also be due to Fabry-Perot (FP) interference in the 
Mica substrate. To rule this out, we repeatedly performed 
Raman spectroscopy on a single 3 QL nanocrystal and sub- 
sequently cleaved the back surface of the Mica to reduce its 
thickness. The resultant spectra all overlapped (not shown), 
implying that FP interference in the Mica substrate can be 
neglected. Nonetheless, for a fixed substrate thickness, the 
interference effects due to multiple reflections in the Bi2Se3 
nanocrystals will change as the crystals are thinned. To check 
the FP effects on the Raman spectra we have performed a cal- 
culation similar to the one shown to work well in graphene and 
Bi2Sr2CaCu2C>8 nanocrystals. 20,21 In Fig. [3] we plot the mea- 
sured and the calculated intensity (FP model) as a function of 
crystal thickness for the main mode (175 cm -1 ). Changes in 
the ratio (4^) due to FP effects are plotted in Fig. [2^ and only 
reveals a very small dependency on crystal thickness, which is 
opposite to what we observe. Therefore, the emergent mode is 
intrinsic to Bi2Se3 nanocrystals and not caused by FP effects. 
However, Fig|3] shows that the change in the overall Raman 
signal with crystal thickness is explained by FP interference 
effects only if the optical constants are modified from the bulk 
values. This modification is not unreasonable given recent 
photoemission experiment a 14 ' 16 . Interestingly, these data also 
show the utility of Raman measurements, as the overall in- 
tensity can be used to determine the thickness of the Bi 2 Se 3 
nanocrystals. 

While the FP interference described above can account for 
the overall intensity-thickness trends in Fig. [T] the origin of 
the 158 cm -1 mode remains unclear. Interestingly, an addi- 
tional mode also appeared in Raman spectra of nanocrystals 
of the isostructural topological insulator Bi2Te3. 8 This mode 
matched perfectly the frequency of an infrared-active mode 
and so it was attributed to the breaking of inversion symmetry. 
Shahil et al suggested that mechanical exfoliation resulted in 
breaks within a QL as well as between them. A similar expla- 
nation may be appropriate for Bi2Se3: a recent IR reflectance 
study of the bulk revealed a mode with the same frequency 
and width 2 ^ (see Fig|2|3). However, in Bi2Te3 this mode ap- 
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FIG. 3. (Color online) Raman intensities as function of crystal thickness at 
175 cm -1 (squares) and calculated Raman intensities due to FP interference 
effects (solid line: optical constants modified from bulk values). 



peared in crystals thinner than 84 nm, whereas it only appears 
in Bi2Se3 nanocrystals thinner than 10 nm. We believe the 
mode may be due to the built-in electric fields at the surface. 
Specifically, the band bending inherent to materials with sur- 
face states will generate asn electric field that will break the 
inversion symmetry. 

We have shown that Raman spectroscopy is an effective 
nanometrology tool for identifying nanocrystals of the topo- 
logical insulator Bi2Se3. This is accomplished by monitor- 
ing i) the overall intensity of the Raman signal and/or ii) the 
strength of the emergent mode at 158 cm^ 1 . The overall 
thickness dependence of the intensity can be accounted for 
by proper modeling of the effect of interference on the Ra- 
man spectra. The origin of the emergent mode remains un- 
clear, although the presented optical conductivity data sug- 
gests that inversion symmetry breaking leads to an IR mode 
becoming Raman active. However, it is interesting to note that 
the emergent mode appears for the same thickness regime of 
Bi2Se3 for which a gap has been theorized to open due to cou- 
pling of the two surfacesJ^i^ In addition, the polar surface 
of Mica may lead to band bending, which breaks inversion 
symmetry. Therefore performing Raman spectroscopy on a 
suspended crystal would provide conclusive determination of 
the influence of the Mica substrate on the Bi2Se3 nanocrys- 
tals. Nonetheless, we have provided a path for fabricating and 
identifying Bi2Se3 nanocrystals through the combination of 
mechanical exfoliation on transparent substrates and the use 
of Raman spectroscopy. This work paves the way for future 
devices and studies of the surface states of topological insula- 
tors. 
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